Thermal spin current and magnetothermopower by Seebeck spin tunneling 

R. Jansen^, A.M. Deac^, H. Saito^ and S. Yuasa-^ 
^ National Institute of Advanced Industrial Science and Technology (AIST), 
Spmtromcs Research Center, Tsukuba, Ibaraki, 305-8568, Japan. 
^ Helmholtz-Zentrum Dresden-Rossendorf, Institute of Ion Beam Physics and Materials Research, 01314 Dresden, Germany. 

The recently observed Seebeck spin tunneling, the thermoelectric analog of spin-polarized tun- 
neling, is described. The fundamental origin is the spin dependence of the Seebeck coefficient of a 
tunnel junction with at least one ferromagnetic electrode. Seebeck spin tunneling creates a thermal 
flow of spin-angular momentum across a tunnel barrier without a charge tunnel current. In ferro- 
magnet /insulator/semiconductor tunnel junctions this can be used to induce a spin accumulation 
in the semiconductor in response to a temperature difference AT between the electrodes. A 
phenomenological framework is presented to describe the thermal spin transport in terms of pa- 
rameters that can be obtained from experiment or theory. Key ingredients are a spin-polarized 
thermoelectric tunnel conductance and a tunnel spin polarization with non-zero energy derivative, 
resulting in different Seebeck tunnel coefficients Sj^ and Sjf for majority and minority spin electrons. 
We evaluate the thermal spin current, the induced spin accumulation and A^/ AT, discuss limiting 
regimes, and compare thermal and electrical flow of spin across a tunnel barrier. A salient feature is 
that the thermally-induced spin accumulation is maximal for smaller tunnel resistance, in contrast 
to the electrically-induced spin accumulation that suffers from the impedance mismatch between 
a ferromagnetic metal and a semiconductor. The thermally-induced spin accumulation produces 
an additional thermovoltage proportional to A/i, which can significantly enhance the conventional 
charge thermopower. Owing to the Hanle effect, the thermopower can also be manipulated with a 
magnetic field, producing a Hanle magnetothermopower. 

I. INTRODUCTION 

The interplay of heat and charge transport is the basis of thermoelectrics, enabling the conversion of heat flow to 
electrical power, and vice- versa. Spintronics concerns the interplay of spin and charge transport and has transformed 
magnetic data storage technology and magnetic field sensing. The connection between these two important fields has 
been established in studies of thermoelectric properties of magnetic nanostructuresi"— . This interplay between heat 
and spin transport, now referred to as spin-caloritronics2iii, has recently gained impetus because the combination of 
thermoelectrics and spintronics offers unique possibilities. On the one hand, it provides a new, spin-based approach 
to thermoelectric power generation and cooling. On the other hand, it provides a thermal route to create and control 
the flow of spin in novel spintronic devices that make functional use of heat and temperature gradients. In addition, 
most spintronic nanodevices involve the application of electrical currents, which create thermal gradients that might 
influence magnetic and spin-related phenomena and thereby device performance and efficiency. This underpins the 
importance of understanding the fundamental interactions between thermal and spin effects. 

A notable recent development is the observation of the spin Seebeck effect by Uchida et al.^^. They found that when 
a ferromagnetic material (permalloy) is subjected to a thermal gradient VT, a spin current is injected into a non- 
magnetic metal (Pt) strip attached to the ferromagnet. This spin current is converted into a voltage proportional to 
VT via the inverse spin Hall effect^^. The name "spin Seebeck effect" suggests it is the spin analogue of the classical 
charge Seebeck effect. The latter can be understood in the following way, noting that the electrical conductance 
depends on the energy of the charge carriers. A thermal gradient across a (non-magnetic) conductor causes a flow of 
electrons with energy above the Fermi energy from the hot to the cold side. Simultaneously, electrons with energy 
below the Fermi energy flow in the opposite direction. There is a net current because the two current components 
do not cancel when the conductivity for electrons above and below the Fermi energy is different. In open circuit 
conditions this results in a voltage between the hot and cold end of the conductor, proportional to S \7T, with S the 
Seebeck coefficient. In a ferromagnetic conductor one expects that the Seebeck coefficient is different for electrons 
with majority and minority spin, as their electronic properties are different. A thermal gradient across a ferromagnet 
would then yield a net flow of spin parallel to the thermal gradient, and produce in a spin voltage (accumulation 
of spin) at the hot and the cold ends of the ferromagnet. Although this was originally suggested to be the cause 
of the observed spin Seebeck effect^'^, the currently accepted interpretation is rather different. It is now considered 
to originate from a non-equilibrium between the magnon distribution in the ferromagnet and the electrons in the 
attached non-magnetic metal, resulting in thermally driven dynamical spin pumping across the interface, without a 
global spin current or spin accumulation in the ferromagnelii^^— . This microscopic mechanism bears no relation with 
the classical charge Seebeck effect. Yet, the spin Seebeck effect is a novel method to convert a thermal gradient into a 
voltage, via the spin, and the phenomenon is generic, i.e., subsequent to the original demonstration for permalloyi^, it 
was also observed in ferromagnetic insulators^, ferromagnetic semiconductorai2i^ and other ferromagnetic metals^i. 
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In a different type of experiment, Slachter et al^ demonstrated for the first time that the Seebeck coefficient of 
a ferromagnet depends on spin. They showed that a thermal gradient in ferromagnetic permahoy induces a spin 
current in the permalloy parallel to the heat flow, and that when the heat flow is directed towards an interface with 
a non-magnetic metal, the spin current crosses the interface and produces a spin accumulation in the non-magnetic 
metal. This thermally-driven spin injection is directly proportional to the difference of the Seebeck coefficient of 
majority and minority spin electrons in the ferromagnet, which was shown to be a fraction of the regular charge 
Seebeck coefficient of the ferromagne t i ^3 ^ 

A distinctly different phenomenon, Seebeck spin tunneling, was observed by Le Breton et ali^. Unlike previous 
work, Seebeck spin tunneling is a pure interface effect that occurs in tunnel junctions with a temperature difference 
AT between the two electrodes, provided that at least one of the electrodes is ferromagnetic. It was demonstrated 
that Seebeck spin tunneling creates a flow of spin angular momentum across a tunnel barrier without a charge tunnel 
current. This thermal spin current was shown to be governed by the variation of the spin polarization of the tunneling 
process with the energy of the tunneling electrons. As will be explained here, Seebeck spin tunneling is directly linked 
to the spin-dependent Seebeck coefficient of a magnetic tunnel contact. This implies that the results of Le Breton et al. 
effectively demonstrated that the Seebeck tunnel coefficient for majority and minority spin is different. In addition, Le 
Breton et al. used Seebeck spin tunneling for thermally-driven spin injection into a semiconductor, i.e., they observed 
that in ferromagnet /insulator/silicon tunnel contacts, the thermal spin current induces a spin accumulation A/i in 
the silicon. 

An interesting analogy exists between electrical and thermal spin transport across a tunnel junction. Seebeck 
spin tunneling is the thermoelectric analogue of spin-polarized tunneling, which refers to the spin dependence of 
the electrical conductance of a magnetic tunnel contact. The latter was clearly demonstrated four decades ago in 
experiment o^^'^^ on ferromagnet/insulator/superconductor junctions, showing that the charge tunnel current between 
a ferromagnet and a non-magnetic counter electrode is spin polarized. In magnetic tunnel junctions comprising two 
ferromagnetic electrodes, spin-polarized tunneling also gives rise to large tunnel magnetoresistance^Zi^, denoting the 
change of the tunnel resistance as a function of the relative orientation (parallel vs. antiparallel) of the magnetization of 
the electrodes. Analogously, Seebeck spin tunneling produces a tunnel magnetothermopower, i.e., a dependence of the 
thermopower of a magnetic tunnel junction on the relative magnetization alignment of the two electrodes. This tunnel 
magnetothermopower (or tunnel magneto-Seebeck effect) has been theoretically predicted^i^^ and recently observed 
by different groups, first in MgO-based tunnel junctions^i^, and subsequently in AI2O3 junctions'^^. Anisotropy of 
the tunnel magnetothermopower was also reported^^. Last but not least, it was predicted that thermal gradients 
give rise to thermal spin-transfer torques in magnetic heterostructures^''^^ and tunnel junctions'^^ and experimental 
evidence for thermal torques has been presented^i^. 

Le Breton et alM described the salient features of Seebeck spin tunneling by numerical evaluation of a free electron 
model. Here we present a phenomenological framework to describe Seebeck spin tunneling in linear response in 
terms of parameters that can be obtained from experiment and analytical or ab-initio theory. Key ingredients are 
a tunnel conductance with a spin polarization that depends on energy, and the spin polarization of the thermally- 
induced electrical transport across the tunnel barrier. An important aim is to establish the connection with a Seebeck 
tunnel coefficient that depends on spin. We evaluate the thermal spin current, the induced spin accumulation and 
A/x/AT, and show that these are proportional to Sj^ — S^f, where Sj^ and Sj^ denote the Seebeck tunnel coefficient for 
majority and minority spin, respectively. We discuss limiting regimes, and point out that the thermally-induced spin 
accumulation increases for smaller tunnel resistance, in contrast to the electrically-induced spin accumulation that 
suffers from the impedance mismatch between a ferromagnetic metal and a semiconductor—"—. We also compare the 
fundamental limits of thermal and electrical spin tunneling. Finally, we demonstrate that the thermally-induced spin 
accumulation produces an additional thermovoltagc proportional to A/i that can significantly enhance the conventional 
charge thermopower. The thermopower can be manipulated with a magnetic field owing to the Hanle effect, producing 
a Hanle magnetothermopower in junctions with only one ferromagnetic electrode. 

II. SEEBECK SPIN TUNNELING 
A. Model 

We consider a tunnel junction with a ferromagnetic electrode and a non-ferromagnetic electrode, typically a semi- 
conductor (Fig. 1) or metal. It is assumed that the tunnel resistance Rtun is much larger than the resistance of the 
electrodes, such that tunneling limits the transport. The ferromagnetic and non-magnetic electrode are characterized 
by so-called spin resistances r/™ and r^, respectively, describing the ratio of the spin accumulation in the material 
and the associated spin current due to spin relaxatior^Sr— . The value of Rtun relative to r{™ and plays an im- 
portant role in the spin transport, as it determines the coupling between the two spin systems. We will assume that 
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Rtun >> i^i"^, which is usually the case when an interface with a (Schottky or oxide) tunnel barrier is formed (the 
resistance-times-area product is 1-10 ^^w? or larger for magnetic tunnel junctions^i^, while r^™ is typically below 
0.01 ^fim? for transition metal ferromagnets^"— ) . Therefore, one does not need to consider the spin accumulation 
and spin-dependent (electrical or thermal) transport within the ferromagnet, including the spin-dependent Seebeck 
effectSS,!^ due to any temperature gradients within the ferromagnet. The spin resistance of the non-magnetic material 
is normally much larger than We do not make any specific assumptions about the value of Rtun relative to r^. We 
thus cover the regime with Rtun > where the spin accumulation in the semiconductor is effectively decoupled from 
the ferromagnet by the tunnel barrier, as well as the regime with Rtun < fs where the coupling to the ferromagnet 
reduces the spin accumulation in the non-magnetic material^^^— . 
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FIG. 1: Energy band diagram of a ferromagnet /insulator/semiconductor tunnel junction. The semiconductor electrode is at 
temperature T„, whereas the ferromagnet is at T/m- A spin accumulation exists in the semiconductor, described by a spin 
splitting A/1 = /i^ — /i"^ of the electrochemical potential. The applied bias voltage V is referenced to the spin-average potentials. 

For a given bias voltage V and temperature difference AT across the barrier, the tunnel current I'^ for each spin 
(ct = '|~, 4,, denoting electrons with magnetic moment, respectively, parallel and anti-parallel to the magnetization of 
the ferromagnet) is: 

/t = G^[v-^^+ L'^AT (1) 
^G^ (v + ^] + L^AT (2) 
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The first term on the right-hand side describes the electrically driven current governed by spin-dependent tunnel con- 
ductances and G^. It incorporates the effect of the shifts of the electrochemical potentials /i'^ in the semiconductor 
due to the presence of the spin accumulation A/i = — (for convenience we have defined A/i in units of volt). 
Note that the spin accumulation typically decays exponentially with distance from the injection contact and that 
A/i denotes the value of the spin accumulation at the interface with the tunnel barrier, as is relevant for tunneling. 
The second term on the right-hand side describes the thermally-induced tunnel current in response to a temperature 
difference, as governed by and L^, which we will refer to as the thermoelectric tunnel conductances (not to be 
confused with the thermal conductance that describes heat flow). We define AT = T„ — T/m, where T„ and Tfm are 
the temperatures of the non-magnetic and ferromagnetic electrode, respectively, and V = Vn — Vfm, where Vn and 
Vfm are the spin-averaged potentials of the non-magnetic and ferromagnetic electrode. 

The total conductances are G = G^ + G^ and L = L'^+L^, and their spin polarizations are Pq — {G^^ — G^) / {G^^ + G^) 
and Pl = (i^ — L^)/{L'^ + L^). The charge tunnel current / and the spin tunnel current Is are then: 

I^I^ + I^ = GV-PgG (^^^+LAT (3) 
L = I^-I^ = PgGV-g(^^^+PlLAT (4) 

The spin current consists of an electrical {Pq G V) and a thermal contribution (P^ L AT), as well as a correction due 
to the A/i that is induced by the (electrical and/or thermal) spin current. The feedback of A/i on the spin/charge 
tunnel current implies that another (independent) relation between A/x and Is is required to obtain a solution. This is 
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provided by the requirement of a steady-state spin accumulation in the non-magnetic material, which implies that the 
spin current Is injected by tunneling is balanced by the spin current due to spin relaxation in the material, integrated 
over its full spatial extent. The spin current associated with spin relaxation is proportional to the spin accumulation. 
We define a spin resistance of the non-magnetic material via: 

Aa* = 2 /, r, (5) 

In our model, is a phenomenological parameter that describes the conversion of the spin current Is that is injected 
by tunneling, into a spin accumulation Afj,. As mentioned before, Afi denotes the value of the spin accumulation 
right at the tunnel interface. This definition of makes no specific assumptions about the spatial profile of the spin 
accumulation in the non-magnetic material, or the formalism used to compute it. If we use the spin-diffusion equation 
and a spin accumulation that decays exponentially with distance from the injection interface with the spin-diffusion 
length Lsd, then the spin resistance of a unit contact area can be expressed as pnLsd, where pn is the resistivity of 
the non-magnetic material^ . This result is frequently used to analyze experimental data, but note that it requires 
introduction of a somewhat unusual factor of two in eqn. ^ . 



B. Spin current and Seebeck spin tunnel coefficient 

Equations (|3]l5|) fully define the system and allow us to obtain the relevant quantities. We first derive a general 
expression for the spin accumulation, valid for electrical (/ ^ 0) and thermal (AT ^ 0) injection, as well as for a 
combination of the two. We shall discuss the case of purely thermal (/ = 0) and purely electrical (AT = 0) driving 
force later on. The solutions for A/i and the spin current in terms of AT and / are: 



2rs 

Rtun + (1 - P'g) rs 
1 

Rtun + (1 - P^) 



- { , " \ [ (^g) Rtun I~{Pl- Pg) 5'oAr] (6) 



[{PG)RtunI -{Pl-Pg)SoAT] (7) 



where Rtun — l/G is the tunnel resistance and 5*0 = —L/G is the charge thermopower (in the absence of a spin 
accumulation; see below). The first term in eqn. ([T]), proportional to /, is the electrical spin current associated with 
the spin-polarized charge current. The second term is the pure spin current due to Seebeck spin tunneling (driven 
by AT) and will be referred to as the Seebeck spin current. The Seebeck spin tunneling coefficient Sst = A/z/AT 
is obtained by setting / = in eqn. and rewriting [Pl — Pg) Sq in terms of spin-dependent Seebeck tunnel 

coefficients = -L^/G^ and S^^ = -I^/G^ for majority and minority spin, respectively. We then obtain an 
important result, namely, A/i/AT is proportional to [Sst — ^st)- 

Sst = = ( D *' , /I d2 A „ 1 (-^It - ^it) (8) 



AT i Rtun + (1 - P<^) r 

Since < Pq < 1, the pre- factor always has a positive sign. The sign of Sgt is thus determined by the difference 
between S'Jj and Sgf Note that the pre- factor tends to zero when the tunnel spin polarization becomes very 
large {Pg ~ 1 or —1), but Sst does not because also G"^ for one of the two spin channels goes to zero. Hence, 
either Sjt or Sst diverges. Taking this into account, one finds that Sst = -\-^{L^ /G'^){rs/ Rtun) for Pg — 1 and 
Sst — —^{L^/G^){fs/Rtun) for Pg — —1- Expressions © and ([7]) apply to situations for which / and AT are fixed. 
This includes recent experiments on Seebeck spin tunneling^i (where / = 0), as well as most experiments on electrical 
spin injection that are performed in constant current mode (and AT = 0). Alternatively, we can express A/i and Is 
in terms of AT and V. However, care has to be taken not to set V = I Rtun, as this is not correct when A/x ^ or 
AT ^ 0, see eqn. ©. 



C. Cliarge thermopower and Hanle magnetothermopower 

The charge thermopower S is obtained from the voltage V\j^q for which / vanishes. From eqn. ^ we obtain: 

V\i=^ _ ^ , (Pg 



^=^ = '5o+(^)5.t (9) 
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This is another important result. In the presence of a spin accumulation, the charge thermopower is not equal to 
Sq = —L/G. There is an additional, previously unidentified, contribution that is proportional to Afj, and thus to 
the Seebeck spin tunnel coefficient Sst- Since Pq can be positive or negative depending on the properties of the 
ferromagnct/insulator interface, and also Sst can have either sign, the additional contribution can enhance or reduce 
the charge thermopower. The enhancement can be significant because Sst can be much larger than Sq, as we will see 
in the discussion section. 

Next we address how the thermally-induced spin accumulation can be detected as a voltage signal. Just as for 
electrically-induced spin accumulation, this can be done via the Hanle effect, which occurs when the spins in the 
semiconductor are subjected to a magnetic field B at a, solid angle 9 with the spin direction^^i^. This causes spm 
precession and consequently a reduction of A/i depending on 9 and on the product of the spin lifetime and the 
Larmor frequency lol = gt^'sB/h, where g is the Lande g-factor, /jg the Bohr magneton and h Planck's constant 
divided by 27r. When the tunnel resistance is sufficiently large {Rtun ^ fs) such that the coupling of the spin 
accumulation to the fcrromagnet can be neglected, spin precession causes a decay of A/i in a Lorentzian fashion^i^: 

AM = 2/,r,^2/,rO |co.2((?) + (10) 

In the absence of any magnetic field there is no spin precession and the spin resistance is . If we keep AT and the 
charge current / constant, apply a magnetic field perpendicular to the spins, and increase B from zero to a value for 
which (jJlTs 3> 1, the spin resistance is gradually reduced from to zero. The A/i then also goes to zero, even if the 
spin current Is that is injected by tunneling is non-zero. This results in the desired voltage change, which is obtained 
from Eqn. dH) and ^ as: 

AVHanle = ^U,=0 " V^r^yyi = (^^^ A/i|„,=o (H) 

An important point is that this expression is valid irrespective of how the spin accumulation is created. In other 
words, also for a thermally-induced spin accumulation, the detected voltage signal AVHanie is given by Pg/'^ times 
A/i, which is the same relation as for electrically-induced spin accumulation'*'^ . 

In the regime where Rtun < fs, the magnitude of A/i is reduced by the coupling of the spins to the ferromagnet, 
but also the functional dependence of A/i on B is modified and eqn. pUj) does not correctly describe the dependence 
on B and 9. However, the maximum and minimum values of A/i for, respectively, = and ujlTs ^ 1 are still 
properly described. Therefore, the amplitude of the spin accumulation can still be correctly obtained from eqn. (|lll) . 
However, extracting the spin lifetime from the line width of the Hanle curve requires a detailed description of the line 
shape taking the interaction with the ferromagnet into account. 

The ability to manipulate the spin accumulation with an external magnetic field (owing to the Hanle effect) also 
means that Sst and hence the charge thermopower S can be controlled by a magnetic field. We define the Hanle 
magnetothermopower Smag as the relative change of the thermopower between its value in zero magnetic field, and 
the value at A/i = that corresponds to ojlTs ^ 1: 

g _ S\uji^=0 — S\i^^r^-S>l _ f -fb \ / 5^stJwL=0 j ^-^^^ 

""^^ S\uiLTs:»i \ 2 / V '^'o / 

The magnetothermopower is mediated by the spin accumulation, has a variation with magnetic field that is governed 
by the Hanle effect, and occurs in tunnel contacts in which only one of the electrodes is ferromagnetic. It is thus 
different from the tunnel magnetothermopower recently observed in a magnetic tunnel junction with two ferromagnetic 
electrodes, which has a magnetic field variation that is controlled by the angle between the magnetizations of the two 
electrodes^Si^. The Hanle magnetothermopower can be very large because Sst can be much larger than 5*0. 

III. DISCUSSION 

A. Origin and definition of Seebeck spin tunneling 

We have seen that Seebeck spin tunneling occurs when the Seebeck coefficient of a magnetic tunnel contact is 
different for majority and minority spin. Le Breton et al,^^ described the salient features of Seebeck spin tunneling by 
numerical evaluation of a free electron model, and showed that Seebeck spin tunneling is determined by the energy 
derivative of the tunnel spin polarization. In this section we will establish the important connection between these 
two notions, and also clarify the definition of Seebeck spin tunneling. 
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First, we note that (Sl^ — S^^) is proportional to {Pl — Pq) and that the spin accumulation induced by Seebeck spin 
tunneling is proportional to {P^ — Pg)- This is easily understood because when / = 0, any thermally-induced current 
(with polarization P^) must be balanced by an equal but opposite electrically-driven current (with polarization Pg)- 
If the tunnel spin polarization does not depend on energy^ all the induced (electrical or thermal) current components 
necessarily have the same spin polarization, and we have P^ — Pg- In that case Sst = 0, the Seebeck spin current 
vanishes for all AT, and a spin current exists only if the charge current is non-zero. To illustrate that Pl = Pg if 
the tunnel spin polarization does not depend on energy, we express G" and in terms of the tunnel transmission 
function D°'{E) integrated over energy E, as^Ml: 



= -- j D%E) {dEf{E, ^i, T)) dE (13) 
L'^ ^-^^ j D^(E)(E- t,) (dEfiE, /X, T)) dE (14) 

where dEf{E,fJ-,T) is the energy derivative of the Fermi-Dirac distribution function f{E,fi,T). When D^{E) and 
D^{E) have the same variation with energy, we can write D'^{E) — x°'Z3(i?), where the coefficients and x'^ do not 
depend on energy. Inserting this in eqns. (fT3)) and (IT4l) wc find Pq = Pl = [x' ~ X^)/{x^ + X^)- Then Pq and Pl 
are independent of E and Pl ~ Pg- Since (Sj^ — S^^) is proportional to {Pl — Pg), we conclude that = S^f if the 
tunnel spin polarization does not depend on energy. This establishes the connection between the energy derivative 
of Pg and a spin-dependent Seebeck coefficient: S'Jj ^ S^f only if the tunnel spin polarization depends on energy. 

There is ample evidence for the energy dependence of the tunnel spin polarization. Indirect evidence comes from 
the decay of tunnel magnetoresistance with bias voltage in magnetic tunnel junctions^^'^^. Direct evidence is provided 
in two reports for transition metal ferromagnets on AI2O3, where the variation of the tunnel spin polarization 
with energy of the tunnel electrons was determined—"^. A significant asymmetry in the decay of the tunnel spin 
polarization with energy below and above the Fermi energy was reported, the decay being much faster above the 
Fermi energy. 

Let us define the term Seebeck spin tunneling more precisely, because one could argue that a thermally-driven spin 
current can exist even if Pl = Pg- While technically correct, in this case the charge current is non-zero, i.e., it is not 
a pure spin current. In fact, for Pl = Pg any thermally- induced spin current is from the spin-polarized charge current 
that arises from the shift of the I — V curve by an amount equal to the charge thermovoltage 6*0 AT. We do not 
consider this to be Seebeck spin tunneling, which is associated with a non-zero energy derivative of the tunnel spin 
polarization (in analogy with the conventional charge Seebeck effect that is related to a non-zero energy derivative 
of the charge conductivity) . Experimentally, one will have a combination of a thermally-driven spin-polarized charge 
current and Seebeck spin tunneling, unless one measures at / = 0, as done in Ref. [2^ . 



B. Comparison of electrical and thermal spin current 



It is instructive to compare the magnitude of the spin current due to electrical and thermal spin tunneling. Besides 
the fundamental interest, this is of course important from a technological point of view. One question is whether 
the creation of a spin current by a temperature difference across the tunnel barrier can be more energy efficient than 
creating a spin current electrically. Another question is whether the heat that is produced by electrical generation of 
a spin current can be re-used to supplement it with a thermal spin current, and how much increase in spin current, 
or reduction in energy consumption, can be obtained in this way. The answer to those questions cannot be given 
in general terms. The efficiency of creating the temperature difference depends crucially on the thermal design of 
the structures. Moreover, whereas it has been known for four decades that the electrical tunnel conductance is spin 
polarized and the polarization has been rather well optimized, the Seebeck spin tunneling has only recently been 
observed and the difference {S^^ — S^f.) is far from optimum. We will therefore only discuss the factors that determine 
the ultimate limits of electrical and thermal spin current. 

We consider the driving term as well as the proportionality factor between Afi and the driving term (see Table 
1). The thermal driving term is 5oAT, which should be compared to the electrical driving term i?f„„ /. For non- 
magnetic metal tunnel junctions Sq has been evaluate d^°'^^ to be in the range of 50 — 100 /iV/K, although it has been 
predicted that it can be enhanced by magnons in ferromagnetic tunnel junctions^i^. The AT for tunnel junctions will, 
in practice, be limited to about 10 K. Hence, ^qAT is of the order of 1 mV, which is to be compared to typical values 
of a few 100 mV for Rtun I- In general, the thermal driving term will thus be smaller than the electrical driving term. 

With respect to the proportionality factor, for electrical spin injection it is limited by Pg, since its absolute value 
cannot be larger than 1 (by definition). However, such a restriction does not exist for the proportionality factor of 
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the thermal spin current, since there is no hmit for the energy derivatives that govern (S'Jj — Sji) and {Pl — Pg)- In 
principle, and can be equal in magnitude but of opposite sign, so that L « and P^ goes to infinity. Physically, 
this corresponds to the situation where the tunnel spin polarizations for states above and below the Fermi energy have 
opposite sign, such that one type of spin is driven from the hot to the cold side of the tunnel contact, and the other 
type of spin is driven from the cold to the hot side. Hence, the proportionality factor for thermal spin accumulation 
can, in principle, be arbitrarily large for suitably engineered materials. This can therefore (more than) compensate 
for the smaller thermal driving term. This suggests that Seebeck spin tunneling can be a viable approach to create a 
spin current, either by itself, or in conjunction with an electrical spin current. 



TABLE I: Comparison of thermal and electrical spin current in a tunnel junction. 



Method 


Type of spin current 


Driving term 


Typical values 


Polarization factor 


Extreme values 


Electrical 
Thermal 


Spin-polarized charge current 
Pure spin current (1=0) 


Rtun I 

SoAT 


~ 100 mV 
~ 0.1 - 1 mV 


Pg 

'-'st '-'st 


±1 

± oo 



C. Magnitude of Seebeck spin tunnel coefficient 

The magnitude of the thermal spin current (and spin accumulation) depends on the value of the polarizations Pl 
and Pg, as well as on the coupling of the spin accumulation to the ferromagnet. An important point is that the 
Seebeck spin tunnel coefficient Sst can be much larger than the regular charge thermopower Sq. To illustrate this, 
the ratio of Sst and Sq is shown as a function of relevant parameters in Fig. 2. For materials with large tunnel spin 
polarization {Pg ~ 1 or —1), a very large Seebeck spin tunneling coefficient is produced if Pl and Pq are unequal, or 
preferably, of opposite sign. This situation would occur for ferromagnet /insulator interfaces that have an almost full 
spin polarization of the tunnel conductance for states at and below the Fermi energy Ep, but a rapidly decaying or 
even opposite spin polarization above Ep, for instance due to the onset of a contribution to the tunneling of a band 
with opposite spin orientation. Since the total thermopower is given by the sum of 5*0 and {Pq/2) Sgt (eqn. ©), the 
thermally-induced spin accumulation in the non-magnetic material can significantly enhance the charge thermopower 
of a tunnel junction. 




FIG. 2; Seebeck spin tunnel coefficient Sst normalized to 5*0, as a function of Rtun/rs, for various values of Pl/Pg and fixed 
Pg = 0.9, and (inset) for fixed Pl/Pg = -2 and Pa varied. Note that Sst = if Pl/Pg = +1. 
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D. Scaling with tunnel resistance 



A noteworthy difference between electrical and thermal creation of a spin accumulation is the scaling with tunnel 
resistance (Fig. 3, and appendix A for explicit expressions for the different regimes). For electrical spin injection, the 
polarization of the injected current {Is /I) is Pg as long as the tunnel resistance is larger than the spin resistance of 
the semiconductor (see appendix A, eqn. (jA4[) ). However, when Rtun ^ rs, the coupling to the ferromagnet starts 
to play a role, and the feedback of A/i on the tunnel transport severely reduces the spin polarization of the tunnel 
current (which is well established^^—). As a result, A/i/I, the spin accumulation per unit injected charge current, is 
constant at large Rtun but decays at small tunnel resistance (Fig. 3, top panel). 

In contrast, the scaling of the thermally-induced spin accumulation is opposite: A/i/AT increases as the tunnel 
resistance is lowered, and reaches a large and constant value when Rtun becomes sma^er than (fig. 3, bottom panel). 
This behavior is consistent with that obtained by numerical evaluation of a free electron model^. It thus appears that 
thermal injection is more efficient at small tunnel resistance, whereas electrical creation of spin accumulation requires 
sufhciently large tunnel resistance to overcome the impedance mismatch. Note that in Fig. 3 we have neglected that 
So decays for thinner tunnel barriers, because the decay is known to be very wea k^"'^^ and does not critically affect 
the main scaling trend. It is also known that for an ultrathin tunnel contact, the thermal (heat) conductance Iq/AT, 
with Iq the heat current, is limited by the interfaces rather than the bulk thermal heat conductance of the tunnel 
barrier material. Therefore, the thermal conductance Iq / AT is expected to be approximately independent of Rtun ■ A 
similar trend would thus result if we would plot A^/Iq instead of A/z/AT. We thus find that the spin accumulation 
per unit charge current is maximum for large tunnel resistance, whereas the spin accumulation per unit heat current 
across the tunnel barrier is maximum at small tunnel resistance. 
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FIG. 3: Scaling of the electrical and thermal spin accumulation with tunnel resistance. Shown are Aji/I for electrical injection 
(top panel) and A^/AT for Seebeck spin tunneling (bottom panel), both as a function of the ratio Rtun/ra of the tunnel 
resistance and the spin resistance of the semiconductor. The results are normalized to the maximum value as indicated (2 Pq rs 
for electrical and Sj^ — S^t for thermal). The dashed lines describe the result when one neglects the feedback of A/i on the spin 
current injected from the ferromagnet. The arrows indicate the reduction due to the feedback. The Pg was set to 0.3. 



IV. SUMMARY 



A phenomenological framework has been presented to describe Seebeck spin tunneling, the thermoelectric analog 
of spin-polarized tunneling. It was established that Seebeck spin tunneling originates from the spin dependence of the 
Seebeck coefficient of a tunnel junction with a ferromagnetic electrode, i.e., Sjt ^ Sjf The connection with a tunnel 
spin polarization Pq that depends on energy was also made. Seebeck spin tunneling creates a thermal flow of spin- 
angular momentum across a tunnel barrier without a charge tunnel current. In ferromagnet/insulator/semiconductor 
tunnel junctions it allows creation of a spin accumulation A/j, in the semiconductor by a temperature difference AT 
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between the electrodes. We expressed the thermal spin current, the induced spin accumulation and A/i/AT in terms 
of the spin-dependent Seebeck coefhcients, tunnel resistance and spin resistance of the non-magnetic electrode. The 
thermally-induced spin accumulation produces an additional thermovoltage proportional to A/z, which can significantly 
enhance the conventional charge thermopower. Because the spin accumulation can be manipulated via the Hanle effect, 
the thermopower depends on a magnetic field, producing a Hanle magnetothermopower in junctions in which only one 
of the electrodes is a ferromagnet. The thermally-induced spin accumulation was shown to be maximum for smaller 
tunnel resistance, in contrast to the electrically-induced spin accumulation that suffers from the impedance mismatch 
between a ferromagnetic metal and a semiconductor. While the efficiency of electrical spin injection is limited by the 
fact that \Pg\ < 1, no such restriction exists for thermal spin current that is determined by the energy derivative 
of Pq, which is unbounded. With suitably engineered materials, Seebeck spin tunneling is thus a viable option for 
efficient creation of spin current. 



Appendix A: Spin current and accumulation in limiting regimes 



1. Thermal spin current 



There are two limiting regimes for Seebeck spin tunneling (/ — 0). When Rtun ^ fs, the induced A/i remains 
relatively small and the feedback of A^ on the tunnel transport is negligible. For this regime one obtains: 

^ « -2 {P, Pa) (-^) So = {l- Ph) [i^) [Si ~ Si,) (Al) 

The spin accumulation decays at larger tunnel resistance since for a tunnel contact 5*0 — —L/G depends only weakly on 
Rtun- This is because all (thermal or electrical) tunnel current components, and hence L and G, decay exponentially 
with tunnel barrier width and height^i^. 
When Rtun/rs < (1 - Pq)^ we have: 

AT ^ [ ^0 = - S^t) (A2) 

G 

In this regime, which corresponds to tunnel contacts with sufficiently low resistance-area product, Sgt and A/i do not 
directly depend on Rtun- 



2. Electrical spin current 



For comparison, the spin accumulation induced by electrical injection of a spin-polarized charge current, 

without a temperature difference across the tunnel barrier, is given by: 



2 rs Rt 
Rtun + {l-Pi;)rs 



Am-^' = ^ 'nl^^. Pg I (A3) 



and the spin current is: 



Rt 

Rtun+ll~P^)rs 



If^^ ^7^^^^^152TTr^G/ (A4) 



We remark that it is customary22ri^ to replace the real tunnel resistance Rtun by (1 — Pq) r*g, introducing r*^ as an 
effective tunnel resistance. The pre-factor then takes a more simple form without the factor (1 — Pq)- As a result, 
it is no longer evident that in order to determine the transition into the regime where the feedback of A/i on the 
tunneling becomes relevant, one has to compare the tunnel resistance to (1 — Pq) Vg. The transition thus depends on 
the value of Pq. We therefore choose to retain the term (1 — Pq) explicitly. 
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